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Internal gravity waves contribute to fluid mixing and energy transport, not only in oceans but also 
in the atmosphere and in astrophysical bodies. An efficient way to transfer energy from large scale 
to smaller scale is the parametric subharmonic instability. We provide here the first experimental 
measurement of the growth rate of this instability. We make careful and quantitative comparisons 
with theoretical predictions for propagating vertical modes in laboratory experiments. 

PACS numbers: 92.05.Bc, 47.35.Bb, 47.55.Hd, 47.20.-k 

Internal gravity waves (IGW) result from the balance of inertia and buoyancy forces in a density stratified 
fluid. Such waves have received a great deal of attention recently because of their relevance and ubiquity in 
different physical situations: they are believed to be of primary importance as they affect ocean mixing and 
energy transport^. Although internal gravity waves do not play the dominant role in the evolution of weather 
and climate, their influence is non- negligible in the dynamics of the atmosphere^. From a fundamental 
point of view, these waves are also particularly intriguing. A striking consequence of stratification is an 
anisotropic dispersion relation relating the frequency to the direction of propagation of the wave and not 
to the wavelength. This property is also encountered for inertial waves (in presence of rotation) or plasma 
waves (in presence of a magnetic field). This has unexpected and interesting consequences in the propagation, 
reflection^ or transmission properties of these waves^. 

Internal waves are known to be inherently unstable due to parametric subharmonic instability (PSI)^. PSI 
is a type of resonant triad interaction where nonlinear terms in the equations of motion allow for efficient 
transfer of energy from large to small length scales where it can be dissipated. The terminology "parametric 
sub-harmonic" is used because, for inviscid fluids, PSI transfers energy from a primary wave to two recipient 
waves of half the frequency. As viscosity effects set in, the frequencies of the recipient waves diverge from half 
the frequency of the primary wave. In previous laboratory experiments, PSI has been qualitatively observed 
by driving low-order standing modes with plungers on the sides of the container^, with an oscillating paddle^i^ 
or relying on the parametric forcing of the tank^. For large amplitude forcing, "irregularities" or "traumata" 
were observed, which led to mixing and overturning. In reS, the critical amplitude of the instability has 
been measured. Quantitative measurements of the growth rate of the instability have never been reported. 

We present here experiments performed with a wave generator, which produces sinusoidal vertical modes 
propagating along a rectangular tank. We measured the growth rate of the instability. This quantity is of 
paramount importance to single out the major mechanism in dissipation processes, a recently highly debated 
issueiii^iii. We first briefly outline theoretical aspects of this instability, after which the experimental 
configuration is described. Then we present our experimental results and compare some of them with 
theoretical predictions. 

Theory Internal waves are characterized by the buoyancy frequency, = V (~^/Po)(dpo/d2:), in which 
g is the acceleration of gravity, po the characteristic fluid density and {dpo/dz) the density gradient in the 
vertical direction z. At large Prandtl number, the 2-D Boussinesq equations of motion can be written as 

where p = p — po is the perturbation density field, i/j the stream function, J the Jacobian operator and v 
the viscosity. Seeking wave solutions with wave number k = (k^m)^ Eq. ([T]) leads to the inviscid linear 
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dispersion relation for frequency uj, 
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For small amplitudes, it can be assumed that that several waves concurrently exist simply as a linear 
superposition. In the case of a resonant triad interaction, where three waves satisfy the spatial resonance 
condition 
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and the temporal resonance condition 
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the nonlinear terms of Eq. ([T]) act as forcing terms transferring energy between the three waves. Each 
wave must satisfy the dispersion relation (|2j). A finite amplitude, large length scale, high frequency wave 
(ko^ujo) can transfer energy to produce two secondary waves of smaller length scales and lower frequencies, 

{ki^uji) and (/c2,cj2)- The instability results from a competition between nonlinear effects and viscous 
dissipation. The growth is exponential if the amplitude of the secondary waves is initially small compared 
to the amplitude of the primary wave^^'^^. In this case, the growth rate is equal to 
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where tpo is the amplitude of the stream function of the primary wave, Ii and I2 are the interaction coefficients 
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and g = 0, 1 or 2 while Ti = fi^'l /2 is the viscous damping factor of the wave i and = /c^ 
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Experimental Configuration A tank, 160 cm long and 17 cm wide, is filled with linearly stratified salt 
water with constant buoyancy frequency N using the standard double bucket method. A monochromatic 
vertical mode-1 wave is generated using the wave generator employed in previous experimental^. The 
generator is composed of 50 plates moving horizontally to impose the horizontal velocity component of a 
mode-1, ie, u{x = 0, t) = —auo^ cos{7tz/H) cos(cjo^), H being the water depth, ujq the excitation frequency 
and a the amplitude of the oscillation of the plates. The motion of the fiuid is captured by the synthetic 
schlieren technique using a dotted image behind the tank—. A camera is used to acquire images of this 
background at 1.875 frames per second. The CIVx algorithmpii computes the cross-correlation between the 
real-time and the t = background images, giving the variation of the horizontal, Px{x^ t) = dx{p{x, t) — 
po{z))^ and vertical, pz{x^z^t) = dz{p{x^ z^t) — po{z))^ density gradients. 
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FIG. 1. (Color online) Snapshot of the horizontal density gradient (in kg/m^) obtained at t = 20T, t = SOT and 
t — 40T with the parameters a — 0.5 cm, cjo = 0.95A/' and N — 0.822 rad- s~^. On the left-hand panel the direction 
of the group velocity, Vg^ is indicated. After 20T, the primary wave has only reached a steady state in the first 20 
cm from the generator (see Fig. [4j). 
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FIG. 2. (Color online) Real part of the Hilbert Transform of px at t = 40T (see Fig. [T]^ right)) presented after filtering 
at cjo = 0.95N (left), cji = 0.38N (center) and 002 = 0.57N (right). Note that the color scale is the same as in Fig.[Tl 



Results Snapshots of an experimental horizontal density gradient field at different times for a particular 
experiment are presented in Fig. [TJ At early times, a pure vertical mode-1 wave can be seen propagating to 
the right away from the wave generator located at x = 0: this is the primary wave. After several frequency 
periods T (typically 30), this wave is destabilized and two secondary waves appear, with different frequencies 
and wave numbers from the primary wave. To see these waves more clearly, the horizontal density gradient 
at later times is filtered at the frequency of the primary wave, CjOq and at the frequencies of the two secondary 
waves uji and UJ2' As described below, the frequencies uji and UJ2 were determined from a power spectrum. 
The result is shown in Fig.[2l Some of the energy of the primary wave has been transferred to both secondary 
waves, leading to a decrease in the amplitude of the primary wave (compare the left part of Fig. [TJleft) and 
Fig. [2]( left)). These two waves have smaller frequency and also smaller wavelength. In agreement with the 
dispersion relation, which links the frequency to the angle of propagation of the wave, the angle of constant 
phase is different for the two wavelengths. For the experiment presented in Fig. [TJ the three measured 
frequencies (co'o, c^i, a;2) are equal to (0.95, 0.38, 0.57)A^, attesting that the temporal resonance condition ^ 
is satisfied. To justify that the spatial resonance condition (|3|) is also satisfied, the components of the three 
wavevectors have to be measured. This is done by extracting the phase of the signal at a given frequency, 
<^o,i,2^ i ^0,1,2^ i ^0,1,2^, using a Hilbert transform^^. At a fixed time and x (respectively z), the phase 
is linear with the position z (resp. x). The component mo, 1,2 (resp. A:o,i,2) is given by the slope of the 
linear fit. Within experimental errors, the wave vectors, represented in Fig. [3l satisfy the theoretical spatial 
resonance condition ([3]). 




FIG. 3. (Color online) Spatial resonance conditions for the experiment presented in Fig.[T] The black line represents 
the theoretical location of the tip of the wave vector ki for a given primary wave vector ko so that the resonance 
conditions Q and ([U and dispersion relations ([2]) are satisfied. The three arrows are the experimental measurement 
of the three wave vectors: the red solid arrow is the primary wave vector ko, the green dashed and blue dotted arrows 
the secondary wave vectors ki and /c2. The square □ represents the most unstable theoretical mode. 



The measured density gradient fields are then analyzed using a time-frequency representation calculated 
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at each spatial point 
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where is a smoothing Hamming window of energy unity Good frequency resolution is provided by a 
large time window h while good time resolution is provided by small time window h. To increase the signal 
to noise ratio, the data is averaged along a vertical line over the water depth. For large ujq/N values, the 
dissipation length is small, so the analysis line is chosen to be close to the generator so that the amplitude 
is large. 

Fig. [Deleft) presents the spectra of the density field for four different excitation amplitudes with ujq = 
0.947V. The spectra are obtained using a time window width equal to 100 T to have good frequency 
resolution. Sx{t,uj) is then averaged over the 10 last periods. Analyzing first the result for the amplitude 
0.5 cm, the picture emphasizes a large peak close to A^, corresponding to the frequency of the mode-1 wave. 
A pair of twin peaks are observed, corresponding to secondary waves of frequencies, uji and CJ2, smaller than 

UJQ. 

The amplitude of each wave is then computed using a time-frequency analysis with a time window width 
equal to 20 T to increase time resolution. The amplitude of the secondary wave of frequency uji is presented 
in Fig. fright). After several forcing periods, a steady state for the primary wave is reached. After a time 
interval, the secondary wave starts to grow and a linear increase of the amplitude on a semilogarithmic plot 
is observed, confirming exponential growth. The value of the growth rate A is measured using a linear fit, 
shown with the dashed lines in Fig. [Upright). The amplitude of the secondary waves eventually saturates. 
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FIG. 4. (Color online) (Left) Spectra of the density gradient field, {Sx{oj))t, for four different excitation amplitudes 
for cjo = 0.94A/' measured at x ~ 12 cm. The amplitudes are respectively 0.25 cm (green), 0.35 (red), 0.5 (black) and 
1 (blue). The value of the amplitude is indicated next to the corresponding curve. The spectra 0.25, 0.35 and 1 are 
multiplied respectively by 0.01, 0.1 and 10 for illustration purposes. (Right) Amplitude of the secondary wave cji, 
*S'a:,i(t), normalized by the amplitude of the primary wave, {Sx,o)t, averaged over time when the steady-state of the 
mode-1 has been reached. Results are similar for the other secondary wave CJ2. The dashed-dotted lines represent 
the amplitude of the primary wave using the same normalization. The dashed lines shows the linear fit, which gives 
the value of the growth rate A. 

Comparing the different curves in Fig. HI one observes that the amplitude has an influence not only on 
the location but also on the height of the peaks of the secondary waves in the spectrum. If the amplitude 
of the primary wave is too small, no peaks are visible and therefore no instability is observed during the 
experiment run time, T^un- This result shows that the growth rate in this particular case has to be smaller 
than 1/run. It may also give an indication of the existence of a threshold in amplitude. As the amplitude 
increases, the distance between the two peaks increases and the instability occurs earlier (after fewer forcing 
periods) and is stronger, i.e. with a larger growth rate which is in agreement with the theoretical growth 
rate (P. 

Experiments were performed using the same stratiflcation and an amplitude of 0.5 cm for frequencies in 
the range of 0.9 < ujq/N < 1. For each experiment, the value of the frequencies of the two secondary waves, 
uji and 6^2, and the growth rate A were measured. Experimental results are presented as a function of the 
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FIG. 5. (Color online) (Left) Values of the frequencies of the secondary waves for the experiments (blue o) and the 
theoretical calculation (red □). The error bars are the half- width of the corresponding peak (vertical error bars are 
smaller than the symbols). The top line of point corresponds to the sum of the frequencies, with corresponding colors 
for theoretical and experimental data. (Right) Values of the growth rate of the secondary waves for the experiments 
(blue o) and the theoretical calculation (red □) . Vertical error bars for the experimental results come from the choice 
of the horizontal location and from the linear fit. These values are smaller than the symbol. T is the period of the 
primary wave. The theoretical values (□) are computed using Eq. ([5]) and the measured amplitude of the primary 
wave. The vertical error bars come from the error made measuring the amplitude of the mode-1. 



frequency of the primary wave, ujq/N^ in Fig. O The sum of the frequencies of the two secondary waves, 
uji + a;2, is equal to the frequency of the primary wave, c^o, within experimental errors, in agreement with 
Eq. ([3]). As CjOq/N increases, the distance between the two secondary frequencies is larger. The measured 
value of the growth rate is presented in Fig.[5fright). The growth rate increases to reach a maximum around 
ujQ = 0.95N and then decreases as ujq gets closer to N. 



To compare quantitatively the experimental results with the theoretical prediction of the growth rate, 
the value of the amplitude of the mode-1 wave has to be precisely known. The theoretical value of the 
amplitude of the streamfunction is equal to aujo/mo. However, the conversion efficiency from the energy of 
the wavemaker to the energy of the mode-1 is less than unity and depends on experimental condition^. 
Moreover, as ujq gets closer to the cut-off frequency, the value of the viscous damping increases^^. 
Consequently, the efficiency is not the same for all primary frequencies cjq, and the amplitude of the primary 
wave has to be measured experimentally to compute the theoretical value of the growth rate. It is important 
to check that the steady-state of the mode-1 wave has been reached. However, the tank being finite in length, 
the measurement has to be performed before the mode-1 wave reflects back into the measurement area. Then, 
using a linear polarization relation, the amplitude ipo of the stream function at this particular frequency 
and wave number is tpo = guoodxpo/ {^k'^pN'^)' The theoretical frequency pair {001^002) of the instability is 
defined as the one that maximizes the growth rate. Without adjustable parameters, the comparison between 
experimental and theoretical results, presented in Fig. O emphasizes a good quantitative agreement. 



Conclusions We have reported the first experimental measurement of the growth rate of parametric 
subharmonic instability in stratified fluids and we have demonstrated this effect in a systematic set of 
laboratory experiments allowing careful comparisons with theoretical predictions. In practice, this heavily 
debated mechanism^^ has implications for many geophysical scenarios. Interestingly, although the generation 
mechanisms of oceanic IGW are quite well understood, the comprehension of the processes by which they 
dissipate is much more open. Consequently, determining the relative importance of parametric subharmonic 
instability, among the four recognized dissipation processes^, is the next step in furthering our understanding 
of how internal waves impact ocean mixing. Quantitative measurements of the subsequent mixing together 
with a fundamental study of wave turbulence would be of high interest. 
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